Theory of extraordinary transmission of light through quasiperiodic arrays of 

subwavelength holes 
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By using a theoretical formalism able to work in both real and k-spaces, the physical origin of 
the phenomenon of extraordinary transmission of light through quasi-periodic arrays of holes is 
revealed. Long-range order present in a quasiperiodic array selects the wavevector(s) of the surface 
electromagnetic mode(s) that allows an efficient transmission of light through subwavelength holes. 
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The phenomenon of extraordinary optical transmis- 
sion (EOT) through periodic two-dimensional (2D) ar- 
rays of subwavelength holes milled in a metallic film [lj 
has sparked a great deal of interest due to both its fun- 
damental implications and its broad range of potential 
applications. Subsequent experimental and theoretical 
works have concentrated on analyzing periodic structures 
0, 0, 13, 0, 0, 0, 11 S E3] • However, very recently, several 
experimental studies showing EOT in quasiperiodic ar- 
rays of holes have been reported 11 , 12, IH, E3] • These 
results suggest that the presence of long-range order in a 
2D hole array is the key ingredient to observe EOT. 

In this Letter we present a complete physical explana- 
tion of the EOT properties found in quasiperiodic distri- 
butions of subwavelength holes. This analysis is based 
on the comparison between the transmission properties 
of finite Penrose lattices of holes with those associated 
with periodic arrays. The picture that emerges from our 
theoretical study is that the physical origin of EOT is 
common for both periodic and quasiperiodic arrays. It 
relies on the excitation of surface electromagnetic (EM) 
modes decorating the metallic interfaces. 

In this paper, our study is focused on analyzing the 
transmission properties of Penrose lattices exhibiting ten- 
fold rotational symmetry, as those studied experimen- 
tally in Ref. 14J. As in the experimental structure, the 
hole radius is chosen to be a — 0.2mm, the thickness of 
the metallic film is h — 0.075mm and the length of the 
rhombus side defining the Penrose tiling, d, is d = 1mm. 
In our simulations, metal is treated as a perfect conduc- 
tor (i.e. with dielectric constant e = — oo), which is an 
excellent approximation in the THz regime. In order to 
calculate the scattering properties and EM field distribu- 
tions, we use a formalism based on a modal expansion 
of the fields at the hole openings 0] , which allows treat- 
ing efficiently large numbers of indentations, arbitrarily 
placed in a metal film. 

Figure 1 shows the three different types of hole ar- 
rangements considered in this work. Left, center and 
right panels correspond to a periodic square lattice, a 
ten-fold Penrose lattice and a random distribution of cir- 
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FIG. 1: (color-on-line) (a-c) Structures considered in this 
work. Square (left), Penrose (center) and random lattices 
(right), (d) Normalized-to-area transmittance (T) spectra for: 
single hole (green line), square array (red line), Penrose lat- 
tice (black line) and a random configuration (blue dots). In 
all four cases, a = 0.2mm, h = 0.075mm and TV = 1506. Inset 
in panel (d) shows the dependence with N for T at resonant 
peaks for the quasiperiodic array, A = 0.83mm (black dots) 
and A = 0.98mm (cyan dots). 



cular holes, respectively. In all three cases, the number of 
holes (N = 1506), their diameter, the film thickness and 
the external radius of the circular array are the same. 
The coordinates in the Penrose lattice were generated 
by the Dual Generalized Method 



15, 16]. The periodic 



structure is a circular portion of a square lattice with 
lattice parameter P = 0.89mm. In the disordered case, 
holes are randomly distributed but without allowing any 
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interhole distance to be smaller than the minimum one 
found in the quasiperiodic case. 

Figure Id depicts the normal incidence transmission 
spectra computed for the three structures, along with 
the transmittance associated with a single hole (green 
line). In all cases, the transmittance is normalized to the 
flux of light impinging on the area occupied by the holes. 
In the spectral range considered, the single hole trans- 
mittance is a smooth decreasing function of the wave- 
length. In the ordered case (red line), the transmittance 
spectrum is also smooth, except close to the resonant 
peak appearing at A = 0.92mm, where the normalized- 
to-area transmittance (T) is about 5 for the geometri- 
cal parameters we are considering. This is the canonical 
EOT peak, appearing at a resonant wavelength slightly 
larger that the lattice parameter. As in the experiments, 
resonant transmission also appears when holes are ar- 
ranged in a Penrose lattice (black curve in Fig. 1). In 
this case, maximum transmission values of about 1.5 are 
obtained at two resonant wavelengths, A = 0.83mm and 
A = 0.98mm. The agreement between theory and experi- 
ment is remarkable (see Fig. 2c in Ref.[14]). On the other 
hand, blue dots in Fig. Id demonstrate that EOT does 
not appear for any distribution of holes: the transmis- 
sion spectrum for the random array does not show any 
resonant feature. This is just a representative example 
of disordered arrays; we have generated several random 
configurations finding always a non-resonant behavior. 

The appearance of EOT can be related to the lattice 
structure in reciprocal space by extending arguments bor- 
rowed from the ordered case to different lattices, as fol- 
lows. Following Ref. 0], EM fields in all space can be 
expressed in terms of the modal amplitudes of the waveg- 
uide modes right at the opening and the exit of the dif- 
ferent holes (E a (R) and E' a (R) , respectively, with R re- 
ferring to the 2D array locations and a running over the 
modes inside the holes). These quantities can, in turn, 
be obtained by solving a coupled system of equations. In 
order to find the link between EOT and the structure 
factor of a given set of holes, S(q) = ^ R exp(— iqR), 
it is convenient to work with the Fourier components 
^a(q) = J2r exp(— iqR)£? Q (R), which satisfy, 

(E„(q) -e n )E n (q)-GXK(q) = J„S(q-ko) 
(SUq)-e„)K(q)-G^„(q) = (1) 

where: 

SW(<0 = E f dk <W S(q - k) £?W(k) (2) 

E n (q) m J 

The expression for the different quantities can be 
straightforwardly obtained from the ones given in Ref. 
Q. External illumination originates I n S(q — ko), ko 
being the in-plane component of the incident wavevec- 
tor. The term e n = —iY n (l + $„)/(l - $„), with 
$„ = exp(2iq n h) (Y n and q n being the admittance and 



the propagation constant of mode n inside the holes, re- 
spectively), is related to the bouncing back and forth of 
EM fields inside the hole. G% = -2iY nV ^/(l - $ n ) 
couples the input and exit sides of the hole. These mag- 
nitudes (e„ and G%) show no dependence on parallel mo- 
mentum, k, as they do not couple modes in different holes 
and they are real quantities for subwavelength holes. The 
terms (q)En} (q) represent the scattering process that 
couples E$ (q) to the continuum Em (k) , the momentum 
difference being provided by the lattice through 5(q— k). 
The amplitude of the process depends on G mn; k: 

G„m ; k = ( 2 \2 E Ywa < "l kcr >-< kcr ' m > ( 3 ) 

where the admittance of the plane wave kcr, Y^ a , is 
g/k z (k) for a p-polarized wave and k z (li)/g for a s- 
polarized one, with g — 2it/\. An important property 
that can be extracted from Eq. (3) is that G mn; k diverges 
whenever a p-polarized diffraction wave goes glancing 
(k z = 0). 
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FIG. 2: (color-on-line) (a) Normalized-to-area transmittance 
versus wavelength for an infinite periodic array (magenta line) 
and several finite square arrays. The geometrical parameters 
are: a = 0.2mm, h = 0.075mm and P = 0.89mm. Inset shows 
the structure factor for the 41x41 case, (b) |Eq — eo| and Go 
versus wavelength for the cases depicted in (a). 

In order to illustrate the mathematics of the formation 
of surface EM modes, let us consider the simpler sys- 
tem of an infinite periodic array of holes. Additionally, 
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we consider normal incidence and assume that only one 
waveguide mode couples directly to external radiation 
(the least evanescent mode, labelled as n—Q). By taking 
advantage of Bloch's theorem [Eq(\c + Gj) = E (\s.) and 
S(k) = J2i <5(k — Gj), being Gi a reciprocal lattice vec- 
tor], Eqs. (1) for k = ko = transform into two simple 
equations for Eq(0) and E' (0): 

(E -e )E (0)-GXE' (0) = I 
(E -eo)E' (0)-G%Eo(0) = (4) 

where So = J^g ^ooiG;- In Figure 2, T (panel a) and 
| So — £ o| (panel b) versus wavelength are depicted for 
an infinite periodic array (magenta line). The geomet- 
rical parameters of this array are the same as the pe- 
riodic one analyzed in Fig.l. As So for G = ±Gi di- 
verges at A = P = 0.89mm, both E (0) and ££(0) 
are zero leading to zero transmission. This is the so- 
called Wood's anomaly [J] or anti-resonance as quoted in 
Ref.Q. The crucial point to realize is that, due to its 
rapid variation close to the divergence, at a wavelength 
slightly larger than the one corresponding to glancing 
angle, |S — eo| = Gq. This leads to a resonant enhance- 
ment of the electric field amplitudes at the interfaces of 
the system [see Eqs. (4)], which can be assigned to the ex- 
citation of a leaky surface EM mode [l7j]. Consequently, 
T presents a maximum at the corresponding wavelength 
(see Fig. 2). Importantly, this resonance appears through 
the coupling to p-polarized modes, closely resembling the 
EM fields of surface plasmons in a real metal. Due to 
that, these modes are usually called spoof surface plas- 
mons emerging when the surface of a perfect conductor 
is periodically corrugated [HI . 

The arguments presented above can be extended to 
the case of finite arrays. Now Bloch's theorem cannot be 
applied and the system of Eqs.(l) must be solved for a 
continuum of states q. However, we have checked that 
for a finite array with a large number of holes q = ko = 
is the dominant transmission channel and the equations 
for Eq(0) and E' (0) could be written like Eqs. (4) with 
S numerically calculated from the knowledge of £*o(k) 
and E' Q (\t) (see Eq.(2)). The results of this approach are 
shown in Figure 2 for the case of square periodic arrays of 
holes (going from 5x5 to 41x41). In this case So presents 
no divergences but resonant features appearing close to 
A = P. The first consequence is that Wood's anomalies 
do not reach zero-value in finite arrays. As for the infi- 
nite case, the cut between |So — eo| and G^ marks the 
location of the transmission peak for large arrays (41x41 
and 31x31). For smaller arrays, there is no cut and the 
transmission peak appears at the wavelength in which 
the distance between |Eo — e o\ and G^ is minimal. 

Once described the periodic case, it is straightforward 
to analyze the case of a quasiperiodic array of holes. In 
panel (b) of Figure 3, the evolution of |So — eo| versus 
wavelength is studied for Penrose lattices with increasing 
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FIG. 3: (color-on-line) (a) Normalized-to-area transmittance 
versus wavelength for several quasiperiodic arrays with dif- 
ferent number of holes, N. The geometrical parameters are: 
a — 0.2mm, h = 0.075mm and d — 1mm. Inset shows the 
structure factor for the N = 1506 case, (b) Both |Eo — eo| 
and Go versus wavelength for the cases depicted in (a). 

number of holes (ranging from N = 106 to N = 1506, 
the case analyzed in Fig.l). |Sq — eo| present maxima 
at wavelengths corresponding to the glancing condition 
for the two main wavevectors of the structure factor (see 
inset of Fig. 3a): b\ (\\ — 0.8mm) and 62 (A2 = 0.94mm). 
Consequently, T shows two minima at these two wave- 
lengths. At slightly larger wavelengths, the distance be- 
tween I So — eo I and G^ is minimal and, correspondingly, 
two transmission peaks appear in the spectrum. There- 
fore, these resonant transmission peaks stem from the 
excitation of surface EM modes at the metallic surfaces, 
much in the same way as in periodic arrays. Notice that, 
however, in the quasiperiodic case, there is no minimum 
wavevector for diffraction (i.e. the structure factor is 
non-zero for wavevectors with moduli smaller than 
see inset of Fig. 3a). This results in diffraction onto ad- 
ditional propagating modes in vacuum (other than the 
zero-order mode), which leads to both smaller resonant 
peaks and less pronounced Wood's anomalies than those 
emerging in the periodic case. 

It is worth analyzing how is the spatial distribution 
of light emerging from the quasiperiodic array. Figure 
4 renders the transmission-per-hole in a Penrose lattice 
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FIG. 4: (color) Transmission per hole (normalized to the single hole transmission) displayed in a color scale for (a) ordered 
case evaluated at A = 0.92mm, (b) Penrose lattice at A = 0.83mm and (c) Penrose lattice at A = 0.98mm. The geometrical 
parameters are the same as in Figure 1. 



of N — 1506 holes at the two resonant wavelengths 
(A = 0.83mm and A = 0.98mm in panels (b) and (c) 
of Fig. 4, respectively). For comparison, panel (a) of Fig. 
4 shows the corresponding distribution for the ordered 
array at the resonant wavelength 0.92mm. In all three 
cases, incident E-field is pointing along the x-direction. 
In the ordered case, due to finite size effects, the maxi- 
mum transmission is located at the center of the structure 
[l9j |. In quasiperiodic arrangements, the transmission- 
per-holc distribution presents a completely different pat- 
tern: it is far from being uniform, showing the appear- 
ance of some holes with high transmission (hot-spots), 
which arc highlighted in the insets of panels (b) and (c). 
Interestingly, in the Penrose lattice, for a given resonant 
wavelength, hot spots show similar local environment. 
However, the existence of hot spots does not imply that 
EOT in quasiperiodic systems is dominated by very local- 
ized resonant configurations of holes. Calculations (not 
shown here) on finite clusters of holes centered at the hot 
spots show an increase of transmittance as a function of 
number of neighbors included in the cluster. This point is 
reinforced by the fact that the resonant peaks observed in 
the transmission spectra of finite Penrose lattices do not 
saturate for small N values (see inset of Fig. Id). Both 
these results are consistent with the interpretation based 
on extended leaky surface EM modes described above. 

In conclusion, by developing a k-space theoretical for- 
malism, we have been able to demonstrate that the reso- 
nant features observed in the transmission spectra of 2D 
Penrose lattices of holes can be explained in terms of the 
formation of surface EM modes at the interfaces of the 
metal film. Furthermore, we have linked the formation 
of these modes to the structure factor of the hole arrays, 
enabling the understanding of the appearance of extraor- 



dinary optical transmission in more general conditions. 
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